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Abstract. We present constraints on the cosmological 
constant Aq and the density parameter Qq from joint con- 
straints from the analyses of gravitational lensing statis- 
tics of the Jodrell Bank-VLA Astrometric Survey (JVAS), 
optical gravitational lens surveys from the literature and 
CMB anisotropies. This is the first time that quantitative 
joint constraints involving lensing statistics and CMB an- 
isotropies have been presented. Within the assumptions 
made, we achieve very tight constraints on both Aq and 
Qq. These assumptions are cold dark matter models, no 
tensor components, no reionisation, CMB temperature 
TcMB = 2.728, number of neutrinos n^ = 3, helium abun- 
dance Yhc = 0.246, spectral index Ug — 1.0, Hubble con- 
stant Ho ~ 68kms~^ Mpc~^, baryonic density fib = 0.05. 
All models were normalised to the COBE data and no 
closed models (fc = -1-1) were computed. Using the CMB 
data alone, the best-fit model has Aq = 0.60 and flo — 0.34 
and at 99% confidence the lower limit on Aq -I- rig is 0.8. 
Including constraints from gravitational lensing statistics 
doesn't change this significantly, although it does change 
the allowed region of parameter space. A universe with 
Aq = is ruled out for any value of flo at better than 
99% confidence using the CMB alone. Combined with con- 
straints from lensing statistics, Aq = is also ruled out at 
better than 99% confidence. 

As the region of parameter space allowed by the CMB 
is, within our assumptions, much smaller than that al- 
lowed by lensing statistics, the main result of combining 
the two is to change the range of parameter space allowed 
by the CMB along its axis of degeneracy. 
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1. Introduction 

Cosmological tests which are sensitive to Aq and fig (the 
normalised cosmological constant and density parameter, 
respectively) can be used to construct likelihood contours 
in the Ag-rio plane. Each test usually has a degeneracy 
such that moving along a curve (which often approximates 
a line) in the Ao-r^o plane leaves the likelihood (almost) 
unchanged. It has long been realised (e.g. Eisenstein et al., 
1999a, b) that the direction of degeneracy of constraints 
from cosmic microwave background anisotropies is roughly 
orthogonal to that of most other tests. Thus, combining 
the constraints from CMB anisotropies with those from 
other cosmological tests can give much tighter constraints 
than either alone. 

Gravitational lensing statistics provide constraints 
which are degenerate such that Aq and flo values are pos- 
itively correlated. This is also the case with cosmological 
tests such as the product of the Hubble constant and the 
age of the universe, the angular size-redshift (standard 
rod) test and the luminosity-redshift (standard candle) 
test. The opposite is the case with constraints derived 
from CMB anisotropies. Thus, it seems natural to conr- 
bine the constraints from Quast & Helbig (1999, hereafter 
Paper I) and Helbig et al. (1999, hereafter Paper II) with 
an analysis of the type performed by Lineweaver (1998, 
hereafter L98), which in itself already provides quite tight 
constraints. 

It is important to note that all three of our analyses 
have fixed all paranreters except Aq and f2o (though in Pa- 
pers I and II an attempt has been made to estimate the 
effect of the uncertainty of the other parameters on the 
derived constraints on Aq and Qq by varying one parame- 
ter by two standard deviations (see Paper I)). Ideally, an 
investigation such as the present one should incorporate 
the uncertainties in all input parameters into the analysis. 
Such a programme is currently under development. 

In this work, we use the most recent CMB data avail- 
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the constraints with the lensing statistics constraints from 
Papers I and II following the procedure outlined in Helbig 
(1999, hereafter Paper III). The plan of this paper is as 
follows. In Sect. 2 we describe the procedure used to cal- 
culate likelihoods in the Ao-ilo plane from CMB data and 
in Sect. 3 we discuss our results and compare them with 
those from Papers I and II. Sect. 4 summarises our con- 
clusions. For a comparison with other recent constraints 
from a variety of cosmological tests, see Paper I. 

Throughout, as in Papers I, II and III, Aq = A/3Hq 
and Qo — SttCpo/SHq refers to the density of matter, 
i.e. not some 'total density' which in our notation would be 
Ao-|-rio (or perhaps including a contribution from pressure 
as well, which we consider to be irrelevant here). The index 
refers to present-day values, since in general these quan- 
tities are time-dependent. (See Paper I for an overview 
of the rest of our notation and general description of the 
gravitational lensing statistics method.) 

2. Constraints from CMB anisotropies 

Gravitational lensing statistics are not very sensitive to 
JIq. However, CMB data can constrain fig more effectively 
(L98). Following Lineweaver et al. (1997, hereafter L97) we 
have calculated probability contours in the Ap-rio plane. 
This method is based on a x^ minimisation: 



X (Ao,i,f^O,; 



E 



(models (Ao,i, i^o,i) - temp^ 



(1) 



where iVcxp is the number of experiments, models is the 
theoretical predicted fluctuation at the multipole range 
covered by the i-th experiment and tempj represents the 
sky fluctuation temperature measured by the i-th CMB 
experiment. Each pair Ao,i,rio,i in the Ao-r^o plane cor- 
responds to a model. We constructed a matrix of mod- 
els and calculated the x^ ^-^d the likelihood associated 
with it, e 2 . The theoretical power spectra were calcu- 
lated with the help of CMBFAST^ (Zaldarriaga, 1998). 
The models depend on a range of parameters. To make 
the test computationally feasible, we fixed all of them ex- 
cept flo and Aq. We consider cold dark matter models, no 
tensor components and no reionisation. No closed mod- 
els (fc == +1) were computed because CMBFAST does 
not yet support this (we are looking forward to the new 
CMBFAST version which will include these models). The 
CMB temperature was set to Tcmb = 2.728, the num- 
ber of neutrinos to n^ = 3 and the helium abundance 
to Ihc = 0.246. The spectral index used was Us = 1.0. 
All models were normalised to the COBE data.^ Finally, 
the Hubble constant and the baryonic density were set 
to Ho = 68kms~iMpc~i and n^ = 0.05. AU these 
values were based on the best literature estimates and 



on the L98 conclusions. Aq and JIq vary in the range 
-0.48 < Ao < 1.48 and 0.06 < ^o < 0.98 with a resolution 
of 0.04. ilo — 0.02 models were not computed since this 
is inconsistent with our value for fib! these models, and 
all outside the examined parameter space, were assigned 
an a priori likelihood of zero. Otherwise, we have used a 
uniform prior. (See Paper I for further discussion.) (Note 
that this is smaller than the range of parameter space 
covered in Papers I and II, but with a finer resolution. 
Initially, wc explored the parameter space as follows: Aq 
and Op vary in the range —3 < Ag < 1 and 0.1 < i^g < 1.5 
with a resolution of 0.1, slightly smaller than the range of 
parameter space covered in Papers I and II but with the 
same resolution. We restricted ourself to the higher reso- 
lution calculations in the smaller area of parameters space 
as outside of this no significant likelihood is present.) 

To compare data and models, the models have to be 
convolved with the window function of the CMB experi- 
ments. The window function delimits the multipole range 
to which the experiment is sensitive. This can be seen as 
a Fourier transform of the experimental beam function 
in the multipole space (White & Srednicki, 1995). On the 
other hand, to compare results from different CMB experi- 
ments, the effect of the window function must be removed. 
This is accomplished by deconvolving both the model and 
the data (see L97 for more details) . The quantities models 
and tcmpj are the deconvolved values. 

We have built up a CMB data compilation that is 
based on L98 and on the web page provided by Tegmark.^ 
We have also added new data from the Tenerife radiome- 
ters and interferometer (Dicker et al., 1999). A list of the 
data used with their references can be found in Table 1. 
The window functions of each of the experiments have 
been also gathered. We have calculated some of them from 
analytical expressions (White & Srednicki, 1995). The rest 
can be found in each of the CMB experiment web pages 
which can be accessed from Tegmark's web page men- 
tioned at the start of this paragraph. 

We do not actually use the COBE points from Teg- 
mark & Hamilton (1997), since the COBE data are used 
internally by CMBFAST. We include them in Table 1 since 
they appear in Fig. 5. On the one hand, CMBFAST nor- 
malises the power spectra to COBE according to the fit- 
ting formula given in Bunn & White (1997). In order to 
take into account the shape of the power spectrum in the 
region of the COBE data, as well as its amplitude, we have 
multiplied the likelihood obtained from our x^ analysis 
(without the COBE points) with the likelihood (again pro- 
vided by CMBFAST using a formula from Bunn & White 
(1997)) of the corresponding power spectra relative to a 
flat power spectrum. 



^ http://www. sns .ias . edu/~matiasz/CMBFAST/cmbf ast .html. 
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Table 1. CMB data used. The window function is centered at I — IcH and drops to half of its central value at Zmin and Zmax, 
except for COBE, where /min and Imax instead indicate the RMS width of the window function. The COBE points from Tegmark 
& Hamilton (1997) are not actually used in our x^ analysis, but are included here since they appear in Fig. 5; see text for details 



Experiment 


ST {f,K) 


+ (mK) 


-(^K) 


I'Uiin 


IcS 


tmax 


Reference 


COBE 1 


8.5 


16.0 


8.5 


2 


2.1 


2.5 


Tegmark & Hamilton (1997) 


COBE 2 


28.0 


7.4 


10.4 


2.5 


3.1 


3.7 


Tegmark & Hamilton (1997) 


COBE 3 


34.0 


5.9 


7.2 


3.4 


4.1 


4.8 


Tegmark & Hamilton (1997) 


COBE 4 


25.1 


5.2 


6.6 


4.7 


5.6 


6.6 


Tegmark & Hamilton (1997) 


COBE 5 


29.4 


3.6 


4.1 


6.8 


8.0 


9.3 


Tegmark & Hamilton (1997) 


COBE 6 


27.7 


3.9 


4.5 


9.7 


10.9 


12.2 


Tegmark & Hamilton (1997) 


COBE 7 


26.1 


4.4 


5.3 


12.8 


14.3 


15.7 


Tegmark & Hamilton (1997) 


COBE 8 


33.0 


4.6 


5.4 


16.6 


19.4 


22.1 


Tegmark & Hamilton (1997) 


FIRS 


29.4 


7.8 


7.7 


3.0 


10 


30.0 


Ganga et al. (1994) 


Tenerife 


32.5 


10.1 


8.5 


13 


20 


31 


Hancock et al. (1997) 


SP 


32.21 


7.44 


4.08 


31 


57 


106 


Gundersen et al. (1995) 


BAM 


55.6 


27.4 


9.8 


28 


74 


97 


Tucker et al. (1997) 


ARGO 


42.01 


6.41 


7.06 


52 


95 


176 


de Bernardis et al. (1994) 


MAX 


43.44 


7.24 


4.94 


78 


145 


263 


Tanaka et al. (1996) 


Python 1 


54.0 


14.0 


12.0 


68 


92 


129 


Piatt et al. (1997) 


Python 2 


58.0 


15.0 


13.0 


119 


177 


243 


Piatt et al. (1997) 


lAC/Bartol 


55.0 


27.0 


22.0 


35 


53 


79 


Femenia et al. (1998) 


MSAMl 


48.42 


11.95 


7.95 


86 


160 


251 


Cheng et al. (1996) 


MSAM2 


59.34 


12.08 


8.23 


173 


263 


383 


Cheng et al. (1996) 


QMAP Fl Ka 


49.0 


6.0 


7.0 


47 


92 


157 


Devlin et al. (1998) 


QMAP Fl Q 


47.0 


8.0 


10.0 


38 


84 


140 


Devlin et al. (1998) 


QMAP F2 Ka 


46.0 


10.0 


12.0 


44 


91 


138 


Herbig et al. (1998) 


QMAP F2 Ka 


63.0 


10.0 


12.0 


81 


145 


209 


Herbig et al. (1998) 


QMAP F2 Q 


56.0 


5.0 


6.0 


58 


125 


192 


Herbig et al. (1998) 


QMAP Fl+2 Ka 


47.0 


6.0 


7.0 


39 


80 


121 


de Oliveira-Costa et al. (1998) 


QMAP Fl+2 Ka 


59.0 


6.0 


7.0 


72 


126 


180 


de Oliveira-Costa et al. (1998) 


QMAP Fl+2 Q 


52.0 


5.0 


5.0 


47 


111 


175 


de Oliveira-Costa et al. (1998) 


Saskatoon 1 


49.0 


8.0 


5.0 


53 


86 


132 


Netterfield et al. (1997) 


Saskatoon 2 


69.0 


7.0 


6.0 


119 


166 


206 


Netterfield et al. (1997) 


Saskatoon 3 


85.0 


10.0 


8.0 


190 


236 


274 


Netterfield et al. (1997) 


Saskatoon 4 


86.0 


12.0 


10.0 


243 


285 


320 


Netterfield et al. (1997) 


Saskatoon 5 


69.0 


19.0 


28.0 


304 


348 


401 


Netterfield et al. (1997) 


CAT 1 


48.44 


7.67 


5.71 


339 


422 


483 


Scott et al. (1996) 


CAT 2 


45.20 


11.02 


8.13 


546 


610 


722 


Scott et al. (1996) 


RING5M 2 


56.0 


8.5 


6.6 


361 


589 


756 


Leitch et al. (1999) 



3. Results 



3.1. CMB results 



Fig. 1 shows the hkclihood e""2- obtained from the x^ 
calculations over our matrix of cosmological models, for 
the CMB data. The contours correspond to 68%, 90%, 
95% and 99% confidence levels, i.e. the area within the 
x% contour level contains x% of the sum of all the like- 
lihood values (one per pixel) in the plot. (See Paper III 
for further discussion of this point, which is important for 
the detailed comparison of different results in the liter- 
ature.) These confidence limits differ from those used in 
L98 in two ways. First, we plot the x% contour as that 
which encloses x% of the integrated probability density 
in the Ao-rig plane, i.e. joint probability contours in Aq 
and flo , whereas those in L98 correspond to the appropri- 



Thus, our contours are naturally larger than those of L98. 
Second, the contours of L98 are actually Ax^ contours, 
which correspond to the appropriate confidence intervals 
if Gaussianity is assumed, whereas ours are 'real' confi- 
dence contours as defined above. (See Papers I and III for 
further discussion.) 

It is well known that the errors quoted for CMB tem- 
perature fluctuation measurements are not Gaussian. Ac- 
tually, for most of the experiments the error bars quoted 
are asymmetric. The observational temperature fluctu- 
ations are usually calculated using maximum likelihood 
techniques, with the value corresponding to the mean and 
the error to the 1-a cutoff. We have taken asymmetric 
error bars into account by using the positive or negative 
error bar according to whether the theoretical value falls 
above or below the experimental value. There is perhaps 



J. F. Macias-Perez et al.: Gravitational lensing statistics with extragalactic surveys. IV 




Fig. 1. The hkelihood function p{D\\o, Qo,^o) based on the 
CMB data in Table 1. (^o represents the parameters other than 
Ao and ilo; ^o, which corresponds to 'nuisance parameters', was 
held constant for all calculations in this paper. See Quast & 
Helbig (1999, Paper I) for definitions and further discussion.) 
All nuisance parameters are assumed to precisely take their 
mean values. The pixel grey level is directly proportional to 
the likelihood ratio, darker pixels reflect higher ratios. In this 
and all subsequent plots, unless noted otherwise, the pixel size 
reflects the resolution of our numerical computations, the con- 
tours mark the boundaries of the minimum 0.68, 0.90, 0.95 and 
0.99 confidence regions for the parameters Ao and ilo and are 
'real contours' in the sense of the discussion in Helbig (1999, 
Paper III) . The diagonal line corresponds to fc = 0; the area to 
the right of this corresponds to spatially closed models which 
were not examined here. The fact that some grey pixels and 
contours are nevertheless in this region is due to finite resolu- 
tion and interpolation, respectively 



of Gaussianity introduces on the constraints on Aq and 
^Q. The L98 confident limits are based on the x^ method, 
which assumes that the likelihood function is a Gaussian. 
Other groups (e.g. Bartlett et al., 1998a, b, 1999a, b) con- 
sider this approximation to be not justified and base their 
calculations on likelihood functions. Unfortunately, likeli- 
hood functions are not provided for all the CMB experi- 
ments and their calculations are based on approximations. 
How to compute this approximation is still poorly under- 
stood. We have opted for a much simpler approach by 
using the likelihood defined above instead of the confi- 
dence limits provided by the x^ statistics. Contours using 
the maximum likelihood approach seem to be larger than 



to lie between the approach of L97 and L98 and that of 
Bartlett et al. (1998b, 1999a).4 

3.2. Review of lens statistics results 

We have repeated the lens statistics calculations of Pa- 
pers I and II at the higher resolution (AAq ~ AOq = 0.04) 
and in the area (-0.48 < Aq < 1.48 and 0.02 < l^o < 0.98) 
used in the CMB calculations. For comparison with Pa- 
pers I and II, these are shown in Fig. 2. 

3.3. Joint constraints 

We follow the procedure outlined in Paper III in comput- 
ing joint constraints. First, to make sure that the cosmo- 
logical tests are not inconsistent with each other, we plot 
the overlap of various contours in Fig. 3. In Fig. 4, we 
present joint constraints formed by the multiplication of 
the corresponding probability density functions. There is 
(at least a small) region of overlap between the CMB con- 
tour and all the lensing statistics contours at 90% confi- 
dence (and of course at 95% and 99% as well), and with the 
JVAS contour (which we consider to be more reliable than 
the optical or joint contours) at 68% confidence. Thus, we 
consider the CMB and lensing constraints to be consistent 
or, at worse, only marginally inconsistent. 

Note that the joint constraints using lensing statistics 
and the CMB data differ only slightly from those using 
the CMB data alone. (It should be remembered that both 
the CMB and the lensing constraints are probably too 
tight since all parameters (except of course Aq and f^o) 
have been fixed for this analysis. With more and bet- 
ter data, both can be expected to improve in the future, 
while improvements in the theoretical models will reduce 
systematic effects. However, since the input parameters 
to the lensing statistics analysis are in many cases better 
understood than those for the CMB analysis, the lens- 
ing statistics constraints are probably more realistic than 
those from the CMB. One should thus not conclude that 
the CMB constraints make lensing statistics superfluous.) 
Nevertheless, the addition of even the current lens statis- 
tics data tightens the upper limit on Aq and the lower limit 
on fio- While the CMB data alone provide perhaps the 
tightest constraints (with the above-mentioned caveats) 
of any cosmological test, they still allow an area of pa- 
rameter space which is ruled out by other cosmological 
tests, among which are lensing statistics. Not only the up- 
per (lower) limit on Aq (^o) is tightened by adding lens 



'Gaussianity' is an issue in at least three different contexts 
with relation to cosmological constraints derived from CMB 
anisotropics. First, the correspondence between Ax^ values or 
fractions of the peak likelihookd and 'real' confidence contours 
as defined above often assumes Gaussianity. Second, not un- 
related, there is the issue of the Gaussianity of the error bars 
of individual experiments. Third is the question whether the 
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Fig. 2. Left: The likehhood function p(D|Ao,no) from optical gravitational lens surveys discussed in Quast & Helbig (1999, 
Paper I), but with a higher resolution and confined to a smaller area of parameter space. (This makes the positions of the 
contours slightly different; see Helbig (1999, Paper III) for a discussion.) Centre; The same but from the analysis of JVAS 
presented in Helbig et al. (1999, Paper II). Right: The same but joint constraints from the two other plots as discussed in 
Paper II 



statistics constraints to those from the CMB, but also the 
best-fit cosmological model shifts to a lower Ag and higher 
flo value. 

The degeneracy in the Aq-Oq plane is such that, in the 
region of non-negligible likelihood, the constraints from 
the CMB alone as well as the joint constraints with lens- 
ing statistics measure approximately Aq + ilo- This region 
is described by the 99% confidence contour, which covers 
the range in Aq of 0.3-0.8 for the CMB alone. In the case 
of joint constraints, the region is shifted to lower Aq values 
as well as slightly increased in size, the exact values de- 
pending on the (combination of) lensing constraints used. 
The corresponding range for Qq is 0.18-0.57, with a simi- 
lar shifting to higher Qq values (and increase in range) for 
the joint constraints. (Of course, the 99% confidence con- 
tour is smaller than the rectangle defined by the ranges of 
Ao and JIq-) 

Taken together, present measurements of cosmological 
parameters definitely rule out the Einstein-de Sitter uni- 
verse (Ao = 0, ilo — 1), "very probably rule out a universe 
without a cosmological constant (Aq = 0) and tentatively 
rule out a flat (Aq + f^o = 1) universe as well.^ A universe 
with Ao w 0.4 and fio ~ 0.3 seems to be consistent with all 
observational data, including measurements of the Hubble 
constant and age of the universe. It should be noted that it 
is really only the CMB data which are indicating a possi- 
bly non-flat universe. Other combinations of cosmological 
tests (e.g. Roos & Harun-or-Rashid, 1999; Turner, 1999, 
and references therein) tend to allow a flat universe within 
the errors. 



"■ This tentative conclusion should be considered with the 
necessary caution. Apart from caveats arising from the limited 
parameter space explored (i.e. all nuisance parameters were 
fixed), the confidence contours cannot be interpreted straight- 



The x^ minimum for the CMB data is obtained for 
ilo = 0.34 and Aq = 0.60. The power spectrum (with the 
data points) for this best-fit model is shown in the solid 
curve in Fig. 5. This is also the best-fit model when the 
CMB constraints are combined with those from JVAS as 
in the centre plot of Fig. 4. The best-fit model for the 
combination of CMB and optical lensing constraints, ei- 
ther with or without the addition of JVAS constraints, 
has 1^0 = 0.38 and Aq = 0.52; this is shown in the dashed 
curve of Fig. 5. 

The comparison values from this work corresponding 
to those in Table 3 of Paper I are presented in Table 2. 

Holding most of the parameters constant is of course 
a weak point of our approach. Obviously, the goal is to 
explore the entire range of parameter space, incorporat- 
ing uncertainties in all parameters, prior information etc. 
This is computationally very expensive. Alternatively, we 
could also think of a multiparameter maximisation me- 
thod which would provide a 'best fit' value for all parame- 
ters, although assigning an error would not be straightfor- 
ward. This might be risky because of possible secondary 
maxima. In fact, our calculations do show a secondary 
maximum, as can be seen by examining the data men- 
tioned in the URL below, although it is too weak to show 
up in the plots. The fact that the secondary maximum 
occurs around Aq = — 1 and f^o — 1 looks suspicious, 
but there is nothing obviously wrong with the behaviour 
of CMBFAST here (M. Zaldarriaga, private communica- 
tion).^ This does appear to be 'real' in the sense that 
it is what the comparison of the data with the CMB- 
FAST power spectra indicate. Of course, there might be 
unknown systematic effects in the former, but as far as we 
can tell, there are no problems with CMBFAST which 
could produce this. On the other hand, it is probably 



Initially, we did find a bug in CMBFAST for fio = 1 and 
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Fig. 3. Tiie 68% (top left), 90% (top right), 95% (bottom left) and 99% (bottom right) confidence contours for each of the data 
sets shown in Figs. 1 and 2. In order of increasing thickness, the curves correspond to Fig. 1 and, from left to right, the plots 
in Fig. 2 

Table 2. Mean values and ranges for assorted confidence levels for the parameter Ao for our likelihoods from this work, for the 
special case f2o = 0.3. This should be compared to Table 3 in Paper I. Only lower limits are given for the case of the CMB 
alone, as the upper limits all lie in the k — +1 area of parameter space, which was not examined. The contours near the fe = 
line are thus not 'real' and should be ignored. Since k — +1 was not examined, no values for the k — case can be extracted, 
as was the case in the corresponding tables in Papers I and II. However, the k — special case has been examined in detail in 
L97. X denotes the fact that there is no intersection of the confidence contour with the Qo = 0.3 line; equal upper and lower 
limits indicate a tangency 
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Fig. 4. The same as Fig. 2 but combined with the CMB constraints from Fig. 1 




Fig. 5. Data points with error bars and the power spectrum 
for our best-fit model based on the CMB data alone or, as the 
values are the same, on the joint constraints from CMB and 
JVAS (solid) and for the combination of the CMB data with 
either the optical surveys discussed in Paper I or with both the 
optical surveys and JVAS (dashed) (again, the values are the 
same) 



not 'real' in the sense that it might disappear when more 
and/or better input data are used or when a more exact 
code than the current version of CMBFAST — especially 
for the relatively poorly explored area of parameter space 
where this secondary maximum occurs — is used. It should 
be kept in mind, however, that there is no a priori reason 
to exclude a secondary maximum. Also, this secondary 
maximum appears in a part of parameter space which is 
ruled out when a few cosmological tests are considered 
simultaneously (see Paper I), so in that sense it is also 
probably not 'real'. 

The CMB data alone do not rule out closed (k = +1) 
models (see also White & Scott, 1996). The probability 
contours are thus compressed near the line that separates 
the open models from the closed ones. This is due to CMB- 



models. Even if these can be ruled out by (some combi- 
nation of) other cosmological tests, it would be useful to 
extend the calculations formed here to include closed mod- 
els, which would allow for an easier interpretation of joint 
cosmological constraints which include those from CMB 
data. 

As mentioned in Papers I-III, to aid comparisons with 
other cosmological tests, the data for the figures shown in 
this paper are available at 

http : //multivac. jb.man. ac.uk: 8000 
/ceres/data_f rom_papers/ 

or 

http : //gladia . astro . rug . nl : 8000 
/ceres/data_f rom_papers/ 

and we urge our colleagues to follow our example. 

4. Conclusions 

We have performed an analysis similar to that of Line- 
weaver (1998), but have used slightly different input data 
and a slightly different statistical technique. We have then 
combined the constraints in the Aq fio derived from the 
CMB with the results of the constraints from gravitational 
lensing statistics presented in Quast & Helbig (1999, Pa- 
per I) and Helbig et al. (1999, Paper II). 

Using the CMB data alone, the best-fit model has 
Aq = 0.6 and flo = 0.34 and at 95% confidence the lower 
limit on Aq + JIq is 0.8. Including constraints from gravi- 
tational lensing statistics doesn't change this significantly, 
although it does change the allowed region of parameter 
space. A universe with Aq = is ruled out for any value of 
ilo at better than 99% confidence using the CMB alone. 
Combined with constraints from lensing statistics, Aq = 
is also ruled out at better than 99% confidence. 

As the region of parameter space allowed by the CMB 
is, within our assumptions, much smaller than that al- 
lowed by lensing statistics, the main result of combining 
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by the CMB along its axis of degeneracy. This axis of de- 
generacy is along a line of constant f2o + Ao, i.e. along a 
line of constant radius of curvature. Indeed, it is close to 
the Oo + '^o = 1 line, which corresponds to a flat universe. 

The CMB analysis favours Aq « 0.60 and n^ ~ 0.34. 
This confirms the long-established result that the Aq = 0, 
fio = 1 (Einstein-de Sitter) model is ruled out by the data 
and supports the recent evidence for a positive cosmolog- 
ical constant (see Papers II and III for a discussion). The 
combination of results discussed in Paper II and Linewea- 
ver (1998) hinted tentatively that flat {k = 0) cosmological 
models are beginning to be ruled out by the data. How- 
ever, to quantify this, it would be interesting to compute 
power spectra for spatially closed {k = +1) cosmologi- 
cal models, since these are not ruled out by current CMB 
data. 

One should keep in mind that these conclusions assume 
fixed values for the nuisance parameters. On the other 
hand, these fixed values correspond to values which are 
(currently) generally accepted and/or values at the global 
X^ minimum in a larger parameter space examined in L98. 

On the one hand, the CMB data provide good evidence 
for a flat universe, since the allowed region of parameter 
space is very small and very near the fc = line. On the 
other hand, as noted in L98, the allowed region is so small 
that a significant departure from fc = is hinted at. (How- 
ever, one should keep in mind that either of these might be 
a consequence of not taking the uncertainties in the other 
parameters fully into account.) Nevertheless, compared to 
other cosmological tests which allow a larger portion of 
the Ao-51o plane, the CMB provides a strong hint that 
the universe is close to being flat. On the other hand, the 
combination of the CMB data and the data from other 
cosmological tests tend to indicate that the best fit might 
actually be achieved for fc < 0, as discussed in Paper I. 
While various tests might individually allow fc = 0, they 
do so for different values of Aq (or, equivalently, Oq which 
of course is 1 — Ao in a flat universe) so that, in combina- 
tion, they provide evidence against a flat universe. 

If the universe does have fc = or is arbitrarily close 
to it, this can never be proven in practice, though our 
confidence in a measurement indicating fc = would be 
inversely proportional to the size of the error bars. On the 
other hand, if the universe is in fact not flat, then this can 
be proven, by reducing the error bars so that the fc = 
case is ruled out. At present, the question of the sign of fc 
or equivalently (assuming a simple topology) the question 
whether the universe is finite or infinite is still an open 
question. On the other hand, there is strong evidence for 
the fact that it will expand forever. 

Assuming the universe is not flat then, since it is rela- 
tively close to being so, to demonstrate that it is not flat 
it will be necessary to reduce the systematic errors in the 
comparison of observations with theory. This can be done 
by incorporating the uncertainties in all input parameters 



mological tests) . Of course, it is also necessary to explore a 
large enough region of parameter space, in all dimensions, 
in sufficient resolution. 

While the joint constraints leave only a small area of 
the Ao-ilo parameter space which fits the observations, it 
should be remembered that neither the CMB nor the lens 
statitistics analyses we have performed incorporates un- 
certainties in the input parameters: all parameters except 
Aq and f2o have been held constant. This is sensible in a 
first-step approach, but of course an analysis of the full 
parameter space should be performed in order to get ro- 
bust constraints on Aq, ^o and the other parameters these 
analyses are sensitive to. Our suspicion is that as a result 
of this the CMB constraints will relax more than the lens- 
ing statistics constraints, so, despite the impression that 
the figures here give that the CMB on its own is power- 
ful enough to constrain the cosmological parameters, one 
should also include gravitational lensing statistics in 'joint 
constraints' analyses (e.g. White, 1998; Eisenstein et al., 
1999a; Tegmark et al., 1998a,b; Eisenstein et al., 1998; 
Webster et al., 1998). This will become more important in 
the future with the completion of large, well-defined grav- 
itational lens surveys such as CLASS (e.g. Myers et al., 
1999). 
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